Scalar field in the Bianchi I: Non commutative classical and Quantum Cosmology 
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Using the ADM formalism in the minisuperspace, we obtain the commutative and noncommu- 
tative exact classical solutions and exact wave function to the Wheeler-DeWitt equation with an 
arbitrary factor ordering, for the anisotropic Bianchi type I cosmological model, coupled to a scalar 
field, cosmological term and barotropic perfect fluid. We introduce noncommutative scale factors, 
I considering that all minisuperspace variables q 1 do not commute, so the symplectic structure was 

modified. In the classical regime, it is shown that the anisotropic parameter /3± nc and the field cj>, 
^ ■ for some value in the \ e [[ cosmological term and noncommutative 6 parameter, present a dynam- 

ical isotropization up to a critical cosmic time i c ; after this time, the effects of isotropization in 
the noncommutative minisuperspace seems to disappear. In the quantum regimen, the probability 
■ density presents a new structure that corresponds to the value of the noncommutativity parameter. 

PACS numbers: 02.40.Gh,04.60.Kz,98.80.Qc 
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The inclusion of matter into the homogeneous cosmologies is usually done with the inclusion of scalar fields, in 



order to study different scenarios, such as inflation, dark matter, dark energy. However, since the early seventies, the 
problem of the appropriate sources of matter to solve an specific scenario [1|,[2|], and its corresponding lagrangian exists. 



In this paper we consider the scalar field and perfect fluid sources as a first approximation in the noncommutative 
quantization regime, and in particular for the Bianchi type I cosmological model, which is the generalization of the 
flat FRW cosmological model. Bianchi cosmological models are the homogeneous but anisotropic generalization of 
the FRW cosmological model, in which the isotropy is not required. In view of the present data that implies a high 
degree isotropy in the Universe, the Bianchi model need an isotropization mechanism to agree with the observation. 
Otherwise they can only model earlier stages of the Universe. In the past different kinds o matter have been used to 
obtain the isotropization. In this paper we explore another possibility to aquire the isotropization in the cosmological 
model: that noncommutativity among the material content (fields and perfect fluid) and the geometric variables of the 
cosmological models produces the expected isotropy. We will explore the effects of noncommutativity in the classical 
and quantum regimes. 

In the last years there has been several attempts to study the possible effects of noncommutativity in the cosmo- 
logical scenario. In particular, in Q], the authors in a cunning way avoided the difficult technicalities of analyzing 
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noncommutative cosmological models, when these are derived from a noncommutative theory of gravity [5!]. Their pro- 
posal is to introduce the effects of noncommutativity in quantum cosmology, by a deformation of the minisupcrspace, 
and it is achieved by means of a Moyal deformation of the WDW equation, similar to the case of the noncommutative 
quantum mechanics 0,0). Some work has been done in this direction, for example in Q the authors study the im- 
plications of noncommutative geometry in minisuperspace variables for an FRW universe with a conformally coupled 
scalar field, using the Bohmian formalism of quantum trajectories also in [lo| | a noncommutative deformation with 
a scalar field in the scalar-tensor theory of gravity was considered. Recently, some work on noncommutative fields 
has appeared, see for example Carmona et al. ll| , where the authors gave the generalization to the noncommutative 
harmonic oscillator and have shown the most general transformation between the coordinates and momenta. In Vakili 
et al. [l2j | a similar procedure is applied to Bianchi I models without matter, and in [13I ] the matter is included. 

In order to obtein a deformed model, there are several ways to proceed, such as [14j |. where the uncentainty relations 
are generalized, which can be of interest, in order to obtain new information in the cosmological scenario and other 
physical quantities. An alternative method that consider deformed special relativity, is considered in [15| . where the 
authors study the effects of the existence of a fundamental length in a cosmological scenario, this is described by what 
is known as the K-deformation [16j |. Both approaches have the same role, to obtain fundamental length, because is 
generally believe that is a natural feature in all theories which have the duty to obtain the quantization of gravity. 

The aim of this paper is to construct a noncommutative scenario for the Bianchi class A cosmological models, 
introducing noncommutative scale factors, considering that all minisuperspace variables q 1 do not commute, so the 
symplectic structure is modified, in order to study the consequences of the noncommutative in the universe, such as 
the isotropization of the anisotropic parameters (3± , so the decay of the field (f> dynamically. The other methods could 
be implemented and carry the analysis if those methods give a better dynamical isotropization of the anisotropic 
parameters. This will be reported in future work. 

The paper is then organized as follows. In Sec. II, we obtain the Wheeler-DeWitt equation including several matter 
contribution, a scalar field, a barotropic perfect fluid and a cosmological term. Also we obtain the corresponding 
classical solutions. In Sec. Ill we present the commutative quantum solutions for the cosmological Bianchi type I, 
considering various scenarios, including the effect to the A cosmological term. In Sec. IV and V, the corresponding 
non commutative classical and quantum solutions are presented, taking that all the minisuperspace variables do not 
commute. Final remarks, are given in Sec. VI. 



II. THE WHEELER-DEWITT EQUATION 



Recalling the canonical formulation of the ADM formalism of the diagonal Bianchi Class A cosmological models. 
The metrics have the form 

ds 2 = — (N 2 - NjNj)dt 2 + e 2n «e 2 *W wV, (I) 

where N and Nj are the lapse and shift functions, respectively, f2(i) is a scalar and /3y(t) a 3x3 diagonal matrix, 
(3ij = diag(/3 + + \/3f3-, (3+ ~ v3/3-, —2/3+), uj 1 are one-forms that characterize each cosmological Bianchi type model, 
that obey duj 1 — |C- k aj j A uj k , Cj k the structure constants of the corresponding invariance group [2|. The metric for 
the Bianchi type I, takes the form 

ds 2 = -N 2 dt 2 + e^e^+^-dx 2 + e^+'^-dy 2 + e 2f V 4 ^dz 2 , (2) 
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The Lagrangian function is given by 



^Total 



-g(R-2A) + 

J-* matter j 



(3) 



we use as a first approximation a perfect fluid and a scalar field as the matter content, in a comoving frame 

Lxotai = \/— g(R — 2A + WnGp + ^d^d^ + V($), 
and using ([2]), we obtain the corresponding Hamiltonian function 



-3S1 



H 



24 



" P Q + p + + p - + 12P * - 48Ae 6 ° + 3847rGM 7 e- 3(7 - 1)n + 12e 6n V(</>) 



(4) 



(5) 



together with barotropic state equation p = jp. The WDW equation for these models is obtained by the usual 
identification, P q n by — i<9 qt < in ©.with q p = (fi, /3 + , /3_, 0). 

Following Hartle and Hawking [181 ] , we consider a semi-general factor ordering which gives 

d 2 d 2 



-3fi 



dfi 2 d(3 2 dp 2 _ 



o2 ^2 £) 

- 12— + Q— - 48Ae 6 ° + 3847rGM 7 e- 3 ^- 1 ) n + 12e 6n V(0) 
o<b l oil 



* = 0. 



(6) 



This equation is not easy to solve in general, therefore we consider several particular cases that yield analytical 
solutions, always considering the scalar field. 

A. Classical Solutions via Hamiltonian Formalism 

To find the commutative equation of motion, we use the classical phase space variables q M = (fi,/3±, (f>), where the 
Poisson algebra for these minisuperspace variables are 



{q",q"} = {q",P q *}=<^, 

and recalling the Hamiltonian equation ([5]), we obtain the classical solutions. 
The classical equations of motion for the phase variables g M and P q » are 



(7) 



n 

P- 
0+ 



{fj,H} = -^c 



-3S2 P 



{/3_,H} = -e- 3 ^P_, 

{/3 +) H} = ie- 3f *P + , 

{^,H} = e- 3O P , 
1 



Pi 



P_ 

p + 

P,/, 



{Po,H} = -e" 3n [-P^ + P 2 _ 
+ [48Ae m + jb ie - 3 ^- 1)n - \2e m V{<p) 
{P_,H} = 0, -> P = ±pi = const. 



12Pd 



{0,H}=, - P + 



±n! = const. 



In order to have an analytic solution, we have consider different cases 



V{4>) = Vq = cte. — > = const 



mi 



(8) 
(9) 
(10) 
(11) 

(12) 

(13) 
(14) 
(15) 

(16) 



4 



we choose the constant in such way for simplicity. In this case, we can associate the potential field with the 
cosmological term, so we have an effective parameter A e // = 48 A — 12Vb for 7 ^ -1, and A e // = 48 A — 3847rGM_i — 
12Vq for the case 7 = — 1. We obtain the relation for Pn 



Pn = ±V a i - A effe 6n + b 7 e- 3 (7-i)" 



(17) 



where af = nf +pf + mf. Employing this relation together with Eq.([5]), for particular epochs of the universe evolution, 
characterized by different values of the 7 parameter, we present the classical solutions in table III A"l 



Case 



-1, X e ff^Q, p_ 1 = M_i 



2 = _K SL 

4 16 



2 
Pi 



7 = 1, A e// <0, p!=M ie - 6n 

q = 



3|A 



Vo I 
4 I' 



a? = n\ + pi + mj + 384ttGMi, 



7 = 1, A e/ / = 0, pi=M ie - eiZ 
&l = nf + pi + + 3847rGMi, 



7 = 0, A e// = 0, po = M e 

b = 384ttGMo, 
af = nf+pl + mf, 



-3£2 



Table \IIA[ Classical Solutions for 7 



Commutative solutions 



Ln 



4af 



leqd* 



8+ = ±|SJ> arctanh 

r t 3 ai 

/3_ = ± I £3- arctanh 

r 3 ai 

, = ± 4 mi t h 



ft = 

/?+ = 
/?- = 



±Ln 



-4af 



2ai 
2ai 
2a 1 



± 



leqol* 

I — arctanh 

3 ai 



= ±| 21 arctanh 

3 ai 

±4= si arctanh 



2ai 



n = 

8- - 



= ±Ln[t"^T], 
= ±Ln [t _ W], 

±Ln[t ^1]. 



2ai 



n = 

B- = 



±Ln 



± 



bpt I ait 

64 4 I' 

±^Ln [ 16a '+ b " t ] 

± jn Ln 16ai+b t 

3ai L t J 

Jp~ Ln r l6ai + bot l 



, 0, and constraints q, a\ and bo. 



In all cases B± are not bounded, excluding isotropization. In the present work we will find out how the noncommu- 
tativity solve this problem. 

III. COMMUTATIVE QUANTUM SOLUTION 

In the following we consider different cases of the potential that allowed us to obtain exact solutions. 

A. V(0) = Vo 



For this case, the potential field is equivalent to a cosmological term, and it is included in the A e g- parameter. Using 
the method of separation of variables, 
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*(ft,/3±,0) = A(n)B(p+)C(P-)V(<j)), in Eq. © we find 

d 2 A „dA 



dn 2 ^dn 



^TrGM^-^^-^^ + Aeffe^ + a 2 ]^ = 0, (18) 

d 2 B 



-v 2 B = 0, (19) 
-/i 2 C = 0, (20) 

r, 2 V = 0, (21) 

where A e g = +12Vo — 48 A is the effective cosmological constant now, and the separation constants are related by 
H 2 = a 2 - v 2 - Urj 2 . The solutions to (p i 120 ]) and (21]) are 



d(3 2 + 
d 2 C 
df3 2 _ 
d 2 V 



B = A x e i!//3 + +B ie - ivf3 +, (22) 
C = A 2 + B 2 e~W- , (23) 
V = A 3 e w + B 3 e" 1 ^, (24) 

where Aj,Bj, j=l,2,3 are integration constants. 

The general solution to equation (|18[) is very complicated to find, so in the following we will solve particular for 
particular cases. 



1. A c ff = 0, 7 7^ 1, and factor ordering Q / 0. 



For this choice of the parameters, the equation (fT5|) is 

d 2 A „ dA 



Q^g + (384^GM 7 e" 3(7 - 1 ) n + a 2 ) A = 0, (25) 



dn 2 

this equation can be written as an ordinary Bessel equation by making the transformation z = zoe -3 ' 7-1 ) , zq 
384ttGM 7 , 



d 2 A / Q \ dA 1 



dz 2 V 9(7-l) 2 y dz 9(7-1) 

whose solution becomes 



^ + WTW^ + a ^ A = ^ (26) 



where p — 3 ( 7 1 _ 1 ) \j ^a 2 — ( 3(7^-1) ^ with ±ip the order of the ordinary Bessel function. 
Then, we find For the complete wave function, we find the equation 

*(n, = z qo J±ip !_*!) ) t A 3 e "' + B 3 e_i,7 1 

(Ai e il//3 + + Bi c-' lvt) +) [A 2 e 1 ^" + B 2 e~ ifif3 -] . (28) 

To obtain the complete wavepacket out of the wavefunction is complicated, that is because all parameters are 
interrelated, so we limit the study to build 3D wavepacket by considering two variables field. This is carry out for the 
case where the parameter 7 take the value -1 (inflation phenomenon). 
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2. A c ff 7^ 0, 7 = 1, factor ordering Q / 0. 
When we introduce the effective cosmological constant A g and a stiff fluid, equation (fT5|) takes the following form 

where Ci = 3847TGM! + a 2 , and making the transformation z = A c ffe en , we find the solutions in terms of generic 
Bessel functions 



A p (z) = z-%Z p (±^y (30) 
A P (Q) = A c-% n Z p (±±V\e 3n y (31) 



with order p = \ \/Q 2 — 4ci and Z p is a generic Bessel function. 

3. A c ff 7^ 0, 7 = — 1, factor ordering Q 7^ 

The other interesting case, corresponds to 7 — — 1 (the equation of state for inflation), proceeding as before, we 
make the transformation z = bc 6 °, where b = 3847rGM_i + A, we arrive to equation 

4H i+ ?)^ + ^ +a2M=o ' < 32) 

the solution is written with the generic Bessel functions 

A(z) = z-3 Z ±ip (±^M , (33) 



with order p — \ \J Act 2 — Q 2 . 

From equations (f2"2l |2"5I |2~4")) we notice that the dependence of the wavefunction on the variables <p and j3± have the 



same structure. In the aim of building a Gaussian weighted wavepacket 20J, [21J, we fix the attention to solutions in 
two of the four variables. For simplicity we choose f2 and the scalar field <f>. We also choose the particular case in 
which the factor ordering Q is zero, 



*(Q,cf>)=Af J°° e- a ("- b )V^J ip QVbe 30 ) dp, 



(34) 



the corresponding plot is presented in figure (HJ, taking the particular values for a=1.5 and b=1.3; and as in the 
Q> fl^t the integral was performed numerically. In this plot the wave function is not peaked around any particular 
values of the variables (fl, (j>). 

B. V(0) = V c- /3< *, stiff fluid 7 = 1 and A = 



The equation © becomes 

d 2 d 2 d 2 5 2 H 



dfl 2 dp 2 d/3 2 d<j> 2 ^dfl 



- 12— + Q— + 3847TGM! + 12V c 



6Q-0</> 



* = (35) 



FIG. 1: 3D plot liP | 2 corresponding to the Gaussian wavepacket equation ([34} in the commuting case, with fj, = v = 



using the ansatz, #(fi,/3 ± ,0) = T(Sl,<j>)B((3 + )C((3-), on Eq. (35|) we find 

1 d 2 T Q cV 12 a 2 ^ 



d 2 B 



T 



v*B = 0, 



d 2 C 



12V e' 



/x 2 C = 0, 



c = 0, 



(36) 
(37) 



d/3 2 ' d/3 2 

where c = a 2 + 3847rGMi, and the constants satisfy the relation /i 2 = a 2 — v 2 . 

There are many changes of variables for equation (136[) that allows to solve it by separation of variables, we choose 
the following one 



x = -6ft + 



y = -20 



(38) 



which is valid for f3 2 ^ 3, then equation (|36|) . using T = X(x)Y(y), is transformed into the following two equations 

d 2 X Q dX 1 



dx 2 2(3 - (3 2 ) dx 12(3 - (3 2 ) 
d 2 Y 2Q/3 2 dY /3 2 m 2 



[l2V e- x + n 2 ] X = 0, 



dy 2 4(3 - !3 2 ) dy 4(3 - (3 2 ) 



Y = 0, 



with n 2 = c — m 2 . Substituting z = e x in (|39[) , we obtain a Bessel differential equation whose solution is 



X(x) = Exp[- 



4(3 - P 2 ) 



x]J± r 



V 



3-/3 2 



-x/2 



(39) 
(40) 

(41) 



with r = 



4(3-/3 2 ) 2 3(3-/3 2 ) 



, the corresponding order. The other equation has the solution 



Y(y) 



A 4 e 4 < 3 ~< 3 "> y + B 4 e 4(3-/3 2 ) y 



(42) 
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with s = y (kS^) + a^Sy- 

For the remaining case, /3 2 = 3 the solution to the differential equation is 

Qn 6bx + ae 6a! ret, , , 

T = F ie -—[e — ~- — e~ — \, (43) 

here k is a non null separation constant, F\ is an arbitrary integration constant, a = c— Q 2 /4, 6 = 12Vo, and the new 
variables a; and y are 

z = - a ,y= - a ■ (44) 



IV. NONCOMMUTATIVE CLASSICAL SOLUTION 

In the commutative model we know that the solutions to Hamilton's equations are the same as in General Relativity. 
Now the natural extension is to consider the noncommutative version of our model, with the idea of noncommutative 
between the four variables (f2 nc , /?± nc , 0), so we apply a deformation of the Poisson algebra. For this we start with 
the usual Hamiltonian ([5]), but the symplectic structure is modify as follow 

{Pn,p±K = {P + ,P-K = o, K,p q 4^ = i, (45) 

{0,0+}* = Mi> {n,j8-}*=i0a, {0-,0+h =i<?3, 
{n,<^ = i04, {0,/3+}, = i0 5 , {<£,/3_}*=i06. (46) 



where the * is the Moyal product [19| , and the Hamiltonian resulting is 



H nc = — — — | -P 2 2 + P^ + p2_ + 12P^ - Ae 6n " c + 12V(0)e 6n » c + b 7 e- 3(7 ~ 1)n " c = 0, (47) 



but the symplectic structure is the one that we know, the commutative one ([7]). There are two formalism to study the 
noncommutative equations of motion, for the first formalism that we exposed have the original variables, but with 
the modified symplectic structure, 

qtc = K,H}„ 

P% c = {P^,H}„ (48) 
and for the second formalism we use the shifted variables but with the original (commutative) symplectic structure 

Qnc = { qnc ! Hnc } ' 

P> = {P£ c ,H nc }, (49) 
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in both approaches we have the same result. Therefore the equations of motion take the form 



^nc 


= {^,H}* = 


- {fl llc ,H nc } 


e -3J2 nc 


Pn, 




12 




P-nc 


= {/?-,H}* 


= {/3-nc, H nc } 


e -3n nc 


-P- + 




1 9 


fi+nc 


= {/?+,H}* 


— {/3+HC7 Hnc} 


e J nc 


-P+ + 


f Pn - -fP„ 


12 


0nc 


= RH}, = 


{</>+nc, H nc } = 


e -3f2 nc p 


^4 
Z 


p - + f p + + t 
z z 


Pq 


= {Pn,H}* 


= |Po,H nc | = 


e ~312„ c 


6Ae 6 ° 


"= + 3(7 - l)b^e~ 


8 


P_ 


- {P-,HK 


= {P_,H nc } = 


= 0, -> 


P = 


= Pi, 


P+ 


= {P+,H}, 


= {P+,H nc } = 


= 0, -> 


P+ = 


= ni. 




= {P0,H}, 


= {P ,H nc } = 


0, 


P^ = 


mi 



if we proceed as in the commutative case we get the solutions showed in the table [TV 



Case 



7 = -1, A e // 7^ 0, P-i 



M_ 



at = nf 



P? 



Kit 
16 : 



7 = 1, X e ff<0, p 1 =M 1 e 



-612 



Pi 



+ m? +384ttGMi, 



3|A- 



Yn.1 
4 I 



-612 



7 = 1, X eff = 0, pi = Mie 
a? = n 2 + p 2 + mf + 3847rGMi 



7 = 0, A e// = 0, po = M e 
b = 384ttGMo, 
a 2 = nf + pf + mf, 



-312 



Noncommutative Solutions 



iLn 



4af 



leqci* 



fPl 



/3 +nc = ± I |i arctanh 
/5_nc = ±1^7 arctanh 
0nc = ±4gff arctanh 





+ f ( 


2ai 








2ai 





fl nc — g Ln 



-4af 



leqd* 



[e qt " 


e 4 / 


2ai 


2 I 



%ni 



4 

,qt 
4~ " 

4 



?mi, 



2\/12 
, a 2 e -qt^ _ ^ 



Pi 



-ni 



2V12 

2\/l2 
-^Pl 



mi, 



mi, 



± 1 211 arctanh 

3 ai 



/3_ nc = iff 1 arctanh 
(6 nc = ±^=— arctanh 



2ai 



L Pl 

+ 



^1 



61 ( 
8 ^ 4 



04 mi 

2 vis' 



aje 
a 2 e 



-qt 



83 „ 
—Pi 



2\/T2 



mi, 



2\/l2 



mi, 



r 0* " 


4 ( e qt 


2ai 


2^4 



+ afe-^ 



fi nc = |Ln [^t] - fpi - + j^mi, 
(3 +nc = ±Ln [t~ft] + f ai - f Pl - 
/3_ nc = ±Ln [t - ft] + %ai + %m - 
= ±Ln [t ^ 



-fni + -# Pl . 



2^12 



mi, 

; m l: 



^ai + ^ni + ^pi. 



fi nc = gLn 



bat_ 

64 



ail 
4 



04 

2 v / l2 



mi, 



4-lLL n ri6ai+bot] , £i A / a 2 , bo£ , ait _ 63 _ 6 
^ai ^ n L t J + 2 V 1 64 + 4 2 P 1 2x/l2 1 



/3+nc — 

4mi Ln 



b t 2 
64 



± 



" 16ai+bot j 



„2 , b t 2 



ai t 



+ — 

4 ' 2 



2 

ni 



2VT2 
fPi- 



mi. 



Table \IVl Noncommutative solutions for, ^ = —1. 1.0. and constraints a, ai and bn 



(50) 

(51) 

(52) 

(53) 

(54) 
(55) 
(56) 
(57) 



We see that the noncommutativity parameters 9\ and 62 could be chosen, with same fine tunning, in order to lower 
the anisotropy at the present time, but in the long term the anisotripy will reappear. The evolution of the isotropy 
(anisotropy) can be appreciate in figure @. 



V. NONCOMMUTATIVE QUANTUM SOLUTION 



As already mentioned, we are looking for the noncommuative deformation of the Bianchi I quantum cosmological 
model. In order to find the noncommutative generalization we need to solve the noncommutative Einstein equation, 
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functions 
1.4 



effective cosmological term =0.1 




time 



functions effective cosmological term =0.01 
2 




time 



20 40 60 80 100 120 



functions effective cosmological term =0.001 
2 




time 



200 400 600 800 1000 1200 



FIG. 2: Plots of /3±nc and <j> field, that appear in the second line in the table IIVI using the values in the parameters ni = 
1, pi = 1, mi = 1 and 8 = 0, 0.05, 0.2, from bottom to top in the figure. The possible isotropization in the parameters f3± nc and 
the decay in the (j> nc field is saw in function of the curvature, but it appears again in this fields after critical cosmic time t c 
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this is a formidable task, due to the highly non linear eharaeter of the theory g fortunately we can circumvent these 
difficulties by following Ref. [4j]. Now we can proceed to the non commutative model, actually we will consider, that 
the minisuperspace variables q 1 = (f2, , /3_ , 0) do not commute, so that the symplectic structure is modified as 
follows 



[q\q 3 ]=i0 ij , [Pi,Pj]=0, [q\Pj]=4 
in particular, we choose the following representation 

[si,p+] = Wi, [n,/3_] = ift,, \p-,p+] = ^3, 

[n,4>] = w 4 , [</>,p + ] = iB s , [<£,/3_] = i0 6 . 



(58) 



(59) 



where the 9\ parameters are a measure of the non commutativity between the minisuperspace variables. The com- 
mutation relation (|58p or (|59[) are not the most general deformations that define a noncommutative field; but this 
relations give us the isotropization of the models that we are working. 

It is well known, that this non-commutativity can be formulated in term of non-commutative minisuperspace 
functions with the Moyal star product * of functions. The commutation relations (|58p can be implemented in terms 
of the commuting coordinates of the standard quantum mechanics ( Bopp shift) and it results in a modification of 
the potential term of the WDW equation [4, 10 1, and one possibility is, for example, 

"nc - tt + yP+ + yP_ + yP , 

0- nc - /9L - |Pn + |P+ - |P„ 

@'\ 9^ „ 

/5+nc -> 0+ - yPfi - yP- - yP , 

a uz, Of. 

0nc - 0-yPo + yP+ + yP- 



(60) 



These transformation are not the most general in order to define noncommutative fields, we can see for example 
Carmona et al. ll|, where the authors give the generalization of the noncommutative harmonic oscillator and have 



shown the most linear transformation between the coordinates and momenta. In Vakili et al. |12j a similar procedure 
is applied to Bianchi I models without matter. 

However, these shifts modify the potential term as follows 

-48Ae 6[n+ ^ p+ + ^ p - + ^ p * ] +384ttGM e -Ki-^+ e -^P + +^-+^4>\ 



\j(n,0±,<j>A) 



_ 12e6 [a + ^ P++ ^P_ + %P,] v / _ ^4 p 



Using the generalized Baker-Campbell-Hausdorff formula 22j , 

e r,(A+B) _ C - V 2 [A,B] eV A 

and the relation between the variables (|59|) . we find the noncommutative WDW (NCWDW) equation 

8 2 d 2 



(61) 



(62) 



dn 2 



dp 2 



+384 7 rGM 7 e- 3(7 - 1)[ °~ 1 T" 9 +-^ 

vU + T 8 a - T a + - T a. 



* = 



(63) 
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To solve the equation we are going to focus our efforts on the same cases that were considered in the commutative 
model. 

Lets start with case 1, A e fi = 0, 7 ^ 1, and factor ordering Q ^ 0. We propose an ansatz of the form 

$(n, 0±,<j>) = e ±i ^e ±il//3 +e ±iM/3 -^(fi), (64) 
and taking in account that e l6 ~^e v * = e 17,e e 7)X , then equation (|63[1 is written as 

d 2 A „dA 



dfi 2 



Q±± + (a? + 3847rGM 7 e- 3 ^- 1 ) n eT§^- 1 H^+9 1 ^ 2 ^ A = ^ (65) 



where the separation constants have the following relation between them, p 2 = a 2 — v 2 — 12rf. This equa- 
tion can be written as an ordinary Bessel equation by making the transformation z = zoe -3 ^ 7-1 -' , with zq = 
3847rGM 7 e =F i^- 1 )( f, « +9ll/+9 ^), 

2 d 2 A / Q \ dA I/2 

z 2 ^r +11+ - — I z— + — — (z 1 

whose solution is 



dz^ 1 ^ ' 9(7-l)V""d* 1 9( 7 -l)2 ^ + « 2 )-4 = 0, (66) 



with order p = 3^-1) ^ 4a 2 — ( 3(7^-1) )" As expected, in the limit 8i — ► we get the commutative solutions, so the 
transition from noncommutative minisuperspace is well defined. 

For case 2, A e ff ^ 0, 7 = 1, and factor ordering Q 7^ 0, the NCWDW equation has the "noncommutative potential" 
U(fi,A cff ,0) = ^Acfte 6 ^-^^-^ --^ 1 ^! +ci). After separating with W = e^e^+e^^- A{fl), we get 



d 2 A . ^d„4 
dtt 2 



Q ^( Acffe6¥9 ' ) + Cl )^°' (68) 



where the constant Ci is the same, as in the commutative case, and f(0j) = =p3(^4?7 + B\v + 02/-O- The solution for 
equation (j6"5|) is written with the generic Bessel function, making the transformation £ = A ge =F3( ' 6 ' 4 ' ,+eiiy+e2 '^e 6n , 

A(0=C^Z P (±l^), (69) 



with order p = \ \/Q 2 — 4ci. 

As in the commutative case, for the simplest factor ordering, Q = 0, simplified solutions are obtained, also the same 
restrictions on A e ff appear. 

For case 3, A e g ^ 0, 7 = —1, and factor ordering Q 7^ 0, the NCWDW equation has the "noncommutative potential" 
U(fi,A e ff,0) = be 6 [ fi ~^ a +~^ 9 -~Tr' 9 *l, with b = 3847rGM_i + A cff . After doing a separation of variables, using 
vj, _ e ±ir;<A e ±ii//3 +e ±i^_^ anc j proceeding as before, we make the transformation £ = be =F3 ^ 4l?+01iy+02At ^e 6n , we arrive 
to the equation 

>d 2 A / Q\ dA 1 



^#H 1+ i^ + ^ +aV=0 ' (70) 
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FIG. 3: 3D plot for (fi, <f>), corresponding to the Gaussian wavepacket equation (172[) in the noncommuting case using 6 = 4and 
fj, = v = 0. 



the solution is written with generic Bessel functions 



A(0 = C^Z±i P [±U 



(71) 



with order p — \ \J 4a 2 — Q 2 . 

The Gaussian weighted wavepacket for the particular factor ordering Q=0 is 



J — i 



e -a(p-b)2 e -^P0J. p 



e 30 e ^p W4 ) df) 



(72) 



and the plot for \^f\ 2 can be see in figures ((5]|3|) taking two values for the noncomutative parameter 8 = 4,0.5 
respectively. Here again we have dropped the P± variables. As in the commutative case, we can to change the 
variable <f> — > (3±, and we obtain similar behaviour. 

The new structure in the probability density can be see in the figures ([5]) and ([6]), in the first along the line of the 
main structure in figures (1) and (3) when = 4, second in the line Q — 1.7 presenting a peak in the variable <f>. 

Finally, for case V(0) = e~' 3 ^, stiff fluid 7 = 1 and null cosmological term, the NCWDW becomes 

' d 2 d 2 d 2 _ d 2 
dti 2 



d(3 2 d(3 2 _ 



12- 



dcj) 2 
* = 0. 



q4t + 3847TGM! + 12V e 6[n ~^ a+ 
oil 



Choosing the wave function 'F = e ±1!y ^+e ±IAI ^ _ JF(r2, <j)) 1 we have 

~v 2 + m 2 + 3847TGM! + i2V e h ( e ')e 6n -^e-^[ /3a " +6 ^] 



d 2 T 3 2 T 
dn 2 deb 2 



'on 



T = 



(73) 



(74) 
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where h(0j) = (30i - §<9 5 )^ + (36> 2 - §<9 6 )/x. 

It is easier to work in the new variables (x,y) given by 
differential equation 



12(3 - /5 2 ) 



d 2 T 

"dx 2 " 



6Q 



dx~ 



12V e h ^)e- x + . 



in place of (f2,^>). Thus we have the corresponding 
4(/? 2 - 3) A 2 T 



T 



2 



dy 2 



2Q 



dF 
dy" 



0, 



(75) 



15 



Psi A 2 




Omega 



FIG. 6: In these plots we present the variation of |^| 2 with respect to Q when 8 — (commutative universe, big plot) and 
8 — 4 (no commutative universe, small plot), along the line of the main structure in figures 1 and 3. 

For solving the equation ((75|) we choose the ansatz JF(x, y) = e( s_ ^ y X(x), where i = > obtaining the following 

equation for the function X, 



d 2 X Q dX 



12V e f ^e- x + n 2 



X = 0, 



(76) 



dx 2 2(3 - (3 2 ) dx 12(3 - (3 2 ) 

where n 2 = c - m 2 , m 2 = 2Q(s - £) + 4(3 - /3 2 )(s - if and f((9j) = h(^) - A (3 - /3 2 ) (s - t) and s is the constant 
that appears in the commutative case. 

This equation is similar to (|39|) . then the solution is similar, and only suffers a shift in the argument of the Bessel 
function 



X(x) = Exp[- 



Q 



"X]J±r 



V 



4(3 -/3 2 ) J \ I/ 3-/3 : 



re f(»0 e -x/2 



(77) 



For the remaining case, f3 2 — 3 the solution to the differential equation is 



(78) 



here k is a non null separation constant, F\ is an arbitrary integration constant, a = c — Q 2 /4, d = 12Voe. m ^*' 
&e m ^\with 



m(0i) = (30i T ^0 6 )f + (30 2 T T ^0 4 «. 

and the variables a; and y arc defined as in the commutative case (equation (|44D ) . 



(79) 



VI. FINAL REMARKS 



In this work by using the equivalence between General Relativity and Hamiltonian formalism, noncommutative 
scenarios are constructed. This was achieved by deforming the minisuperspace for the Bianchi type I cosmological 
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model coupled to barotropic perfect fluid, field <j) and cosmological term, in the gauge N = 1, as we can see the 
solution f2 nc is the commutative solution plus a function on 0j, independent of time. We have also analyzed the 
f3± nc noncommutative solutions, in some ranges on the parameter (9.; and effective cosmological constant, occurs a 
dynamical isotropization, i.e., (3± nc are bounded in the present time, noting that different evolution with respect to the 
commutative (3±, also the field <j> decays dynamically to a constant, but after a critical time t c the evolution is the same 
in both scenarios. In the quantum stadium, in all but one of the cases considered the influence of the noncommutativity 
is encoded in a change of the argument of Bessel functions and the change is not qualitatively very remarkable ( see 
figures (|1I3I6[) . The singular case is the one with 1 = 3 and here the influence of the noncommutativity is in the 
argument of an exponential in the form 6dx + ae 6x . Only the constant d depends on the noncommutivity parameters 
and therefore we have that for x large and negative they will dominate and for large positive x the noncommutative 
effects disappear. From the definition of x we see that for finite <f> the noncommutative regime corresponds to Q — > — oo 
and at late times commutativity is recovered. It would be of interest to improve the work with the two approaches 
cited in the references [Lll and obtain the possibles differences in the dynamics of our universe, using anisotropic 
cosmological models; this will be reported in future work. 
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